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We find the Holographic Renormalization Group equations for the 
holographic duals of generic gravity theories coupled to form fields 
and spin-^ fermions. Using Hamilton-Jacobi theory we discuss the 
structure of Ward identities, anomalies, and the recursive equations 
for determining the divergent terms of the generating functional. In 
particular, the Ward identity associated to diffeomorphism invariance 
contains an anomalous contribution that, however, can be solved ei- 
ther by a suitable counter term or by imposing a condition on the 
boundary fields. Consistency conditions for the existence of the dual 
arise, if one requires that a Callan-Symanzik type equation follows 
from the Hamiltonian constraint. Under mild assumptions we are able 
to find a class of solutions to the constraint equations. The structure of 
the fermionic phase space and the constraints are treated extensively 
for any dimension and signature. 
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1 Introduction 

In theories with diffeomorphism invariance the Poincare generators are gauged, 
and the concept of a point has no invariant meaning. It becomes therefore all 
but impossible to define local observables. The only known way around this 
problem is provided by Holography |l|, |2[ , which means (in the strong form) 
that a theory with diffeomorphism invariance should be describable in terms 
of a dual local quantum field theory defined on the boundary of space-time. 
In this article we investigate the conditions that a diffeomorphism invariant 
theory of interacting fermions and form fields has to satisfy in order that its 
holographic dual, in the sense defined presently, exists. 

The AdS/CFT conjecture P| has provided a very concrete example of 
the Holographic Principle, asserting that the physical content of a theory 
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of gravity in d + 1 dimensions should be encoded in a quantum conformal 
field theory living on the boundary of its anti de Sitter space-time. The 
correspondence has been made precise in 0, [| where it was explained how to 
compute field theory observables (Green functions) in terms of AdS integrals. 
The central object in this correspondence is the generator of connected Green 
functions W 

e -W[g,J] _ j p $ e -S[ff,*]-J©[#] 

where J is, on the QFT side, a set of external sources coupled to (composite) 
operators O and g^ is a background metric, which allows to compute the 
stress-energy tensor. The conjecture then states that 

W[g,J] = S SVGRA [g,J] (2) 

where the on-shell supergravity action is evaluated on an AdS solution, and it 
is a functional of initial values J of the fields appearing in it. The quantum 
field theory content is therefore encoded in a classical field theory for its 
sources. 

This correspondence has been later generalized to non-conformal theories 
II 0- H, H- |K| 0, 0- There the classical radial equations of motion on the 
supergravity side were interpreted as the renormalization group (RG) flow. 
In this way supergravity solutions provide information about RG trajectories 
connecting different CFT's appearing at the fixed points of the flowQ. In this 
bulk/boundary correspondence the central object is still W, which is now a 
functional of the sources at a given mass scale. 

In QFT W obeys a set of Ward identities, which ensure that classical 
symmetries are preserved at the quantum level and, consequently, constrain 
the counter terms that can be used in renormalization. For instance, the 
Ward identity of (broken) scale invariance is the Callan-Symanzik equation, 
which is a particular form of RG equation. Much of the information con- 
tained in this structure can be traced back to divergences. Regularization of 
these infinities introduces a scale, and the requirement that the renormalized 
quantities are independent of the cutoff produces the RG equation. In this 
article, we will discuss how the holographic renormalization group (HRG) 

3 This mechanism was discussed in |l3| in the framework of Polyakov's proposal for 
holographic noncritical strings |l4|, [l5| . 
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manages, generically, to produce these features, but also how, and why, it 
may sometimes seem to fail if the theory includes interacting fermions. 



As Holography relates theories in d and d + 1 dimensions it is useful 
to resort to Hamiltonian formalism, where the transverse "time" coordinate 
plays a distinguished role. As we are particularly interested in the on-shell 
action it is rather natural, as proposed in , to use Hamilton-Jacobi theory. 
There are also some conceptual reasons for doing so: Most importantly, on 
the QFT side, one eventually needs a formulation in terms of first order 
equations, as they can be interpreted as RG flow equations. Also, on the 
gravity side, it is natural to expect the Hamiltonian formulation to arise 
JT7J , as one considers classical gravity as a WKB limit of quantum gravity 



(String/M-Theory) . 

In the Hamiltonian formalism one chooses a coordinate t, with tangent 
vector field eo, to parameterize the evolution of the initial value hypersurface, 
the "boundary" . It is sufficient to characterize this hypersurface by giving its 
normal direction n with norm n-n = ±1. Depending on this sign, denoted ry, 
one is considering either time-like or space-like boundaries. This will produce 
sign changes with respect to usual Minkowskian Hamiltonian theory, and we 
will most efficiently keep track of them by leaving rj unspecified. Otherwise 
we follow the standard Hamiltonian reduction. Writing the full metric 

ds 2 = (r]N 2 + NiN^dt 2 + 2N i dtdx i + g ij dx i dx j (3) 

the canonical (ADM) gravitational action coupled to matter has the form 

S = j dt jp q - J d d x NH± + JV*Wi + 9 a £a} (4) 

where pq is a shorthand for the kinetic term of the dynamical degrees of 
freedom, N and iV* are the lapse and shift functions, which together with 
Q A are the Lagrange multipliers for the constraints^ 

H ± ~H t ~g A ~0. (5) 

These constraints generate, in the bulk theory, orthogonal deformations and 
diffeomorphisms on the initial value hypersurface and possible additional 
local symmetries, such as gauge, local Lorentz or supersymmetries. 

4 These are first class constraints. We will consider also second class constraints. They 
appear in the presence of fermionic fields that enter the action linearly. 
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In the Hamilton-Jacobi theory one makes a canonical transformation such 
that the new phase space coordinates are constants of motion. In gravita- 
tional theories the generating functional F of the canonical transformation 
must satisfy simultaneously the constraints 

n±(q^) = (6) 

Wifo^O = (7) 
5F 

G A (q,-^) = 0. (8) 

The Hamilton principal function F[q] is actually just the classical action 
evaluated at some given time t for fixed boundary values q(t). The momenta 
can be calculated from 

This gives us the action; if we want to calculate the full equations of motion, 
we have to consider also 

q = {q,H} DB , (10) 

where {-,-} D b denote the Dirac brackets. The holographic correspondence 
allows us to give a dual field theory interpretation to the equations above. 

Halving the degree of the equations of motion by using Hamilton-Jacobi 
equations for F, in spite of giving a rather appealing form for application 
to holography, boils down to having a high degree of arbitrariness in F. 
Equations (^), (0), and (§) have in fact a huge number of solutions, many 
of them unacceptable on physical grounds. This was discussedf] in [jnj , for 
a scalar field coupled to gravity. One should keep in mind also the fact that 
not every gravity solution corresponds to deformations of a CFT by adding 
an operator to the fixed point Lagrangian, as was pointed out in pU| , |2T|. 



Rather, some of them may correspond to an altogether different vacuum 
of the same theory, where an operator has acquired a nontrivial vacuum 
expectation valueQ. In general, only a subset of solutions will admit a direct 

5 There are some comments in this direction in |Q, as well. 

6 For a clear exposition of this point see for instance [ p2[ . This was illustrated explicitly 
in fll9(l in connection to the Hamilton-Jacobi formalism. 
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physical interpretation pfl . To ensure a physically acceptable result one 
specifies the asymptotic (near the AdS boundary) behaviour of the solutions. 
If one wishes to solve the equations for a coupled system, with a dynamical 
metric or interacting sources, in a perturbative fashion, the back-reaction 
will involve the flow equations, which determine the scaling behaviour of the 
fields, modifying then the general ansatz for F — even for its local part, 
which is the only one computable in general. 

Finally, let us briefly comment on how the QFT divergences arise on the 
gravity side. It is well known from explicit calculations that the Einstein- 
Hilbert action (with the cosmological constant and the Gibbons-Hawking 
term) is divergent when evaluated on asymptotically AdS solutions. For 
instance in the coordinate system 

ds 2 = i(dt 2 + 4dxW) (11) 

the determinant gives y/g = t~ d \/g^, and Ry = t 2 R^ } so that divergences 
arise in the limit t — > (asymptotically AdS boundary) for terms containing 
up to [^rp] Ricci's. The existence of such solutions in this context is a physical 
requirement: It means that the RG has fixed points. In general, including 
additional terms in the bulk action will introduce further divergences. 

In this article we shall consider the fermions and form fields with some 
apriorously fixed anomalous dimensions at the UV conformal point, and even- 
tually consider how general the results are. This is tantamount to restricting 
the analysis to a particular subset of operators used to flow from the UV 
CFT. In even boundary dimension there can be additional logarithmically 
divergent terms that, in the Hamilton- Jacob i context, can arise only from 
nonlocal contributions. One should isolate the terms in F that are divergent 
in the CFT limit: These terms will give rise to the beta functions, while the 
rest will give the renormalized generating functional. 

It turns out that - in analogy with what is achieved by QFT Ward iden- 
tities - these local divergent terms, which we will call Sdw, are fixed by 
the constraints in the classical theory, and can indeed be determined recur- 
sively order by order according to their degree of divergence. This recursive 
procedure can be expressed as "descent equations" and is computationally 
equivalent to the counter term generating algorithm proposed in p4| , |25fl in 
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the context of pure AdS gravityQ Performing this calculation one sees that 
(0) and (§) will produce Ward identities (diffeomorphism and gauge symme- 
tries), while (|6|) provides a recursive equation for S^iv, formally very similar 
to the perturbative expansion of the master equation, that in addition will 
give rise to a Callan-Symanzik-type equation for W. 

One should naturally try to chose a renormalization procedure that pre- 
serves the symmetries. A failure to fulfill this requirement would give rise 
to an anomaly. Anomalies in the boundary theory may arise from bulk con- 
tributions to the Hamiltonian constraints. This may happen because the 
Hamiltonian "time" -slicing might not preserve local bulk symmetries after 
gauge fixing. An example is provided by the holographic chiral anomalies 
discussed in [El, and in [E7[] in the Hamilton-Jacobi framework. 



After the work of |16[ there has been a number of papers discussing dif- 
ferent aspects of the HRG using the Hamilton-Jacobi formalism for gravity. 



Most of them, however, elaborate gravity coupled to scalar fields [28, 29 



DJ [H], |32|, [0| |34], |35| , where the situation is by now well understood. In 
this article we shall extend this analysis to essentially more general situations. 
We consider a general setup in which spin-| fermions and form fields interact 
with gravity. In particular the treatment of spinors is rather subtle. These 
fields are linear in momenta, couple derivatively to the metric and are in gen- 
eral endowed with a complex structure. It was already noticed in previous 
literature [^, [37], [H], on the AdS/CFT correspondence that fermion fields 
split naturally in momentum and coordinate parts, and that the latter part 
are the relevant sources for the boundary operators. The split boils down, in 
general, to imposing some Lorentz-invariant condition on the fermions, such 
as chirality or reality conditions. In the following we systematize these obser- 
vations in the context of Holography, using the Hamilton-Jacobi approach. 
The Ward identities resulting in this case from the gravitational constraints 
generalize the structures that arise from the scalar and Yang-Mills sectors, 
and display some intriguing novel features, as well. 

The plan of the paper is as follows. In Section 2 we use the pure gravity 
system to illustrate how to implement and solve the descent equations that 
follow from the 7i± constraint. In the process we will present a novel way 



7 In that context form fields were considered in |26|. 
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of calculating the holographic Weyl anomaly of We make then general 
comments on how Ward identities and holographic anomalies arise. In Sec- 
tion 3 we briefly review the structure of the holographic Callan-Symanzik 
equation and notice that there exists a distinguished coordinate system in 
which the flow equations take the familiar form involving beta functions. Sec- 
tion 4 is devoted to explaining the Hamiltonian treatment of spin-| fermions 
and form fields. Sections |5|, ^, and [7| contain our main results. In Section |5| we 
present the bracket structure induced by the TC± constraint, while in Section 
H we discuss the other Ward identities. In Section [7] we solve, under certain 
(mild) assumptions, the bracket equations. We find agreement with results 
concerning the pure gravity sector and some interesting constraints in the 
form field and fermion sectors. The last section contains conclusions and fu- 
ture perspectives. In the appendices we have relegated some useful formulae 
and notation, and an extensive analysis of the fermionic phase space. 



2 Renormalization and Anomalies 

In order to illustrate the procedure, we consider in this section pure gravity 
with a cosmological constant. We start by showing how to get a solution 
for local terms of the generating functional F. As a by-product we shall get 
an alternative derivation of the holographic Weyl anomaly computed in |3~9|. 
We then set the ground to discuss Ward identities and anomalies. 
Equation (§) in this case reads 

-VK 2 ( 1 \ 5F 5F p: ( 1 „ . \ 

9ik9ji ~ -, — ~9ij9ki -= — t — = V9 ~ A (12) 



V9 V d - 1 J dgij bg H \k z 

and, following the notation introduced in [}40|| , we rewrite it as[] 



(F,F)=C + C 2 . (13) 

We then expand S<n v according to the degree of divergence of each local term. 
For the case of pure gravity this is equivalent to a derivative expansion, as is 
easily seen by a change of coordinates x« — > t~ x Xi. 

The local terms that diverge at the UV fixed point are those that contain 
up to d — 1 derivatives; logarithmically divergent terms may arise in the non- 
local part of W. However, the most divergent ones come always multiplying 

8 A similar approach has recently appeared in |32j. 
9 The subscript counts the number of derivatives. 
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a scale invariant term. This is typical of the structure that arises in effective 
actions that produce a conformal anomaly j|T] . 



Let us set F = S<n v + W in (|T2| ) and expand order by order in t. In global 
scalings the metric behaves as — > t~ 2 gij. Given that 

(,-) - t d (;-) (14) 

Sn _ r d+n S n (15) 
W -> iy-21ogt5 d , (16) 

and (So, S^) = 0, the following descent equations have to be satisfied 

£(D = (So, S ) 

C 2 = 2(So,S2) 
= 2(S ,S 4 ) + (S 2 ,S 2 ) (17) 

-2(S ,W) = 2(S 2 ,S d _ 2 ) + --- + 2(Sa_ 2 ,S i+2 ) + (S^Sa) . 

2 2^ 2 2 

The zeroth order equation fixes the relationship between the cosmological 
constants in the bulk and on the boundary, see Eq. (|7T|). The n'th order 
equation is an equation for S n in terms of S 2 , . . . , S n - 2 . In particular, the 
last equation^ then gives the trace of the bare stress-energy tensor. This is 
because the lowest order term So acts through the bracket (■, ■) essentially as 
a scale transformation. The trace anomaly at the conformal points, where 
the beta functions vanish and couplings go to their fixed point values J*, is 
given by 

<T> = Mlzll lira -L (S 0) MO- (18) 

In the UV limit this expression remains finite and reproduces the holographic 
Weyl anomaly of . The results for d = 2 and d = 4 areQ 



(T) = —jrR (19) 

K 4 A 

(T) = h- (3R^R i3 ~R 2 ) . (20) 



10 In even boundary dimension. In odd dimensions it fixes also higher order terms, and 
nonlocal logarithmic terms do not arise, as is well known. 



"With T] = 1, see equations ([n]) and (f7|). 
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The only Ward identity that one has to check in this case is diffeo- 
morphism invariance, Eq. (0). It is trivially satisfied for the finite (bare) 
stress-energy tensor, because the counter terms are generally covariant in 
the boundary metric 

5W 

v '%r- Vi7X=0 - (21) 

If there are matter fields, apart from the additional Ward identities related to 
other local symmetries, even the diffeomorphism identities will be modified. 
In fact, in the presence of nonzero vacuum expectation values the CFT stress- 
energy tensor is not conserved, but obeys rather, for one scalar source, say, 

V j (Ty) = - (O) Vi J . (22) 

This should not be considered, however, as a diffeomorphism anomaly, but 
rather a feature of the spontaneous symmetry breaking induced by the nonzero 
vacuum expectation values of O. This is in accord with the non-conservation 
of the Brown- York quasi-local stress-energy tensor, as was previously 



noticed in pi |. The signature of a gravitational anomaly is, instead, stress- 
energy nonconservation with no vev's turned on, analogously to what hap- 
pens if there is a holographic chiral anomaly [[], T7\. For obtaining an explicit 



anomaly in pure gravity it would therefore be necessary to consider the ap- 
propriate Chern-Simons terms. Such a term would be, for instance, in the 
(2, 0) theory in 6 dimensions the reduction on S A of the M-theory CS term 

S cs = J d 7 x Tr(cj A R A R A R) . (23) 

We shall see in Section |] that in the presence of fermion fields this pattern 
will be somewhat modified. 



3 The Holographic RG Equations 

In this section we review how the holographic RG equations arise. Here, we 
do not restrict to pure gravity theories, but assume that the brackets in (p~3| ) 
have been modified to accommodate variations with respect to matter fields. 
Equation (|T3") reads 

(Sdi VJ Sdiv) + 2(S div , W) + (W, W) = £ . (24) 
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The second term on the LHS is a linear operator acting on W. (W, W) is a 
quadratic correction to the linearized RG, and the left-over piece gives addi- 
tional non-homogeneous terms, whose leading contribution was shown in the 
previous section to give the Weyl anomaly. In fact, because of cancellations 
achieved in the descent equations fll7f) , fl24j) can be rewritten 

2(S div ,W) + (W,W) = 0(l) . (25) 

In order to get a Callan-Symanzik type equation, one can take the follow- 



ing steps ||16|| : Calculate n variations of ( p5|) w.r.t. the sources and drop the 
contact terms containing more then one delta function, which would not con- 
tribute for operators evaluated at different points. Then letting the sources 
go to constant (in Xi) values, the metric assuming the form, 

9ij = n(t)Vij > ( 26 ) 
and eventually integrating once over dx d one gets expressions of the form 

(upJL + Pj^ (O h O h ) + ^d h (3 J < (OjO, h ) = 0(t d \og 2 t) , (27) 
whereQ 

Pi = (^GiJ S 4f) (28) 

U * ■ (29) 

VV^ °9ij J j=j{t), 9 =n 

Notice that these equations are taken at a finite cutoff t, i.e. for bare quanti- 
ties. Would one allow the couplings to be space-time dependent the vanishing 
of the beta functions would correspond to the sources obeying d-dimensional 
Einstein equations that would follow from S^iv This could give rise in prin- 
ciple to boundary backgrounds that are not Minkowski, or its conformal 
completion. 

Considering finally the Hamiltonian flow equations (0) it turns out that 
in the appropriate gauge the transverse coordinate plays the role of a param- 
eter for the scale transformation induced by the metric. In general it is not 



12 



We drop a numerical coefficient. 
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always true that a gravity solution has the interpretation of a RG flow: It 
depends in fact on the leading behaviour of the bulk fields near the boundary 
of AdS |^0|, |21| . In the Hamilton-Jacobi context this means that not every 



solution Sdiv of the constraints gives rise to physical flows [|T^]. However, for 



those that are physical, the flow equations read, after fixing the gauge Ni = 

ft ~ NU n (30) 
jj ~ NPl (31) 

The scale transformation depends parametrically on the cutoff t. Let us now 
use (0) to express ( p7|) in terms of the variation of t, and choose the gauge 
N = +j. This yields 

n 

+ J2d Ii f3 J t(O h ---0 Ji ---0 In ) = 0(t d log 2 *), (32) 
and (PTf) can be written as 



Pi = t^J z , (33) 

which is consistent with the beta functions defined before and with interpret- 
ing ( |32"D as a bare RG equation. Furthermore the right hand side represents 
higher order terms that, in this spirit, are logarithmic corrections to scaling. 
This choice of gauge differs from the Fefferman-Graham [J5 by a sign: The 
latter corresponds in fact to N — — |. 

As is the case in QFT, the definition of the beta functions is not unique. 
They are fixed unambiguously near the fixed points, where their behaviour 



is universal, but their extrapolation at intermediate RG steps is not |35 
The definition given in [jlfj differs from that given above by the ratio of the 
rate of scale change (j3CD . In fact, in that reference N = 1 and all quantities 
depend on \i. The fixed points are not affected, because the zeros of the beta 
function cannot be modified by dividing by U. Also, it is well know that 
the bare RG equations and the renormalized ones have the same physical 
content. 
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4 Fermions and Form Fields in Hamiltonian 
Theory 

Let us now turn to a generic (d + l)-dimensional theory with spinors, form 
fields and local diffeomorphism invariance. We will consider the most gen- 
eral two- derivative action (neglecting Chern-Simons terms) with quadratic 
fermion couplings consistent with gauge symmetry, namely S = Si + Sn + 
Sin, where 

Si = — J d d+1 x^/g (R - 2r/V • V n n + 2r/V ■ (n tr K) - k 2 a) (34) 
Sn = fd d+1 x^g(^-F A F A + F A J A ) (35) 



\2X 2 

Sin = ^Jd d+1 x^(^Mfl^-(ptlj)Mtlj + 2^Z A T A ^). (36) 

The fields appearing in these expressions are the following: The field strength 
F = dA is an Abelian p-form and couples to the fermions through J A = 
%[)(,T A ip. The capital Latin letters refer to a multi-index of pertinent rank; 
summations include division by the factorial of the rank. There can be arbi- 
trarily many fermion flavours, but we always suppress the index that would 
distinguish them. This action encaptures, and generalizes, many interesting 
features of the effective superstring actions. For instance, F could be thought 
of as a Ramond-Ramond field. 

The nondynamical couplings (, Z A and M mix fermion flavours, and are 
not assumed space-time constants, unless explicitly indicated. They satisfy 
suitable hermiticity conditions so that the action is always real. They can be 
thought of as the Yukawa couplings to higgsed scalar fields. As Z A is not a 
dynamical field, we have assumed Z 0A = 0. The underlying space-time has a 
boundary to which we can associate a normal vector field n and the extrinsic 
curvature K. This vector field has constant norm n • n = r\ = ±1, and is 
hence either temporal or spatial. 

In order to go over to the Hamiltonian formalism, we need to choose a 
particular "time coordinate", or "evolution parameter" whose tangent field 
we call cq. The bulk vielbeine e - decompose into boundary vielbeine L i -, 
the vector n, and the nondynamical degrees of freedom and N l . The 
connection on the boundary is obtained (see Appendix [A]) from that in the 
bulk by shifting the Christoffel symbols in such a way that the n becomes 
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a covariantly constant vector field on the boundary, and that the boundary 
basis is covariantly constant in the normal direction 



Vi n = (37) 
n-ViCj = . (38) 

Given the time direction eo we can now reduce all tensor fields in terms 
of fields defined on the boundary. Isolating the kinetic terms in the action, 
i.e. terms involving derivatives along eo , the remainder is by definition the 
total Hamiltonian 

S = J dt (p q - H) (39) 

where 

pq = Jd d x (> d t L i& + E A d t A A + X 9tV ~ M x) (40) 

H = J d d x (NH± + N'Hi + Ajg A + e^J^j (41) 

We have added here, by hand, the constraint that guarantees the freedom to 
choose the flat basis for vielbeine freely, i.e. the generator of local Lorentz 
transformations in the bulk, J7"— . For notation, that is standard, see Ap- 
pendix 0. 

The physical phase space can now be easily read off from the above refor- 
mulation. It consists of the canonical pairs (p l a ,Lj), (E^,A^), (x ,^); an d 

{ip a , Xb)- The fields A$ , N and N* are Lagrange multipliers that correspond 
to the first class constraints Q A , Ti±_ and Hi, respectively. Due to fermions, 
that have a first order action principle, there are also second class constraints 

X = ^r"Mf (42) 

In order to solve these constraints, we have to split the fermion phase space 
in some Lorentz invariant way into two parts, one of which we treat as the 
configuration space, and hence boundary fields, and the other of which are 
then the momenta, to be solved as functionals of the boundary fields in HJ 
theory. This is done in Appendix |C| by imposing a chirality condition; here 
we shall only point out some salient features. 
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Most of the calculations must be done case by case, but the emerging 
structures are similar. Having used the condition of choice to put the kinetic 
term in such a form (as in formulae ( |131| ) and ( |137| )) that the coordinates 
and the momenta can be read off, there are always some left-over pieces: 

1) There turns out to be a total time-derivative term d t G. This term 



should be simply subtracted from the action, as argued in |L7], |38 
Generally speaking the reason for this procedure is that the generating 
function arises also on the bulk side, eventually, from a path integral, 
which is going to be defined more fundamentally in Hamiltonian lan- 
guage. It has also been pointed out |5B| that, as first order actions 
vanish on classical solutions, we would otherwise get a trivial result as 
far as the fermions are concerned. 

2) There will also arise terms that involve a time derivative of the metric. 
This will naturally change the gravitational momentum, but in a way 
that is easily kept track of. After some algebra it turns out to be 
sufficient to simply shift by 

vr iJ — ► n 11 - l -g l3 G . (43) 

5 The Holographic Callan— Symanzik Equa- 
tion 

Having split the fermion phase space we are now in the position to write 
down the Hamilton-Jacobi equations for the full system. We will perform 
this analysis for Weyl fermions and assume that M is a space-time constant 
matrix. 

It turns out that, provided there are no marginal operators Z present in 
the bulk and that the rank of the tensor field p is odd, the Hamilton-Jacobi 
equation originating from TC± does indeed take the generally expected form 

(F,F)=C. (44) 

Where now 

(F,F) = (F,F) g + (F,F) A + (F,F)^, (45) 
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and the RHS of (M) is 

1 _ . 1 



C = ^ / 9^- 2 R-A + ^F A F A + F A ^T A ¥ , 

+ l -^MIf) V - l -{If)^)M V + ^Z A Y A ^j . (46) 
It is useful to define the following operators^ 

V = 2% (47) 
TP = /— ^ (48) 

ogij 2 

V A = -L + ^CM- 1 ^ + ^M" 1 ^ (49) 

0/1^ o<£> 0(y9 

Had we also considered p even, the last equation would have been different: 
Then the operator T A would have contained terms with either no or two 
derivatives w.r.t. the fermion fields. Now the brackets can be written easily^ 

(F,H) g = ^( gil g jk --J— gij g kl )(T«F)(T kl H) (50) 
(F,H) A = 0,(T A F)(T A H) (51) 

The reason for the fact that also the fermionic momenta give rise to a bracket 
is easily seen in the case for Weyl fermions: The chiralities of the coordinates 
and the momenta are such that if we insert any operator even in Clifford 
matrices between them, 7tO even (p, the result is nontrivial. Similarly, the non- 
trivial results for odd operators arise from insertions between either two 
coordinates, (fOo^tp, or two momenta, nO^ii. 

This structure changes slightly if the bulk mass terms, or couplings to 
external form fields, Z, are included. If Z is even, such as a mass term, there 
will be an additional term 

6 Z F + (F,F) = C (53) 



13 Derivatives act from the left, but not on fields included in the same operator. 
14 Here we have set Z to zero. See below. 



15 



where 



) 



(54) 



This addition preserves the flow equation form of the final Callan-Symanzik 
equations, however. Also, if Z is odd or M is not constant, there will be an 
additional quadratic piece in the fermion momenta, and the basic form of 
the brackets will be, eventually, unchanged. 

If we are forced to give the initial data in terms of Weyl fermions, as we 
are assuming in this article, the requirement that the bulk action produce a 
QFT generating functional that obeys the Callan-Symanzik equation implies 
that bulk theories where fermions are coupled to dynamical even rank form 
fields do not posses a simple holographic dual. 

Including higher order interactions could be a problem, because they 
would introduce higher powers of momenta, which would spoil the basic form 
of the brackets. However, we have seen above an encouraging rearrangement 
of terms, where the structure of the theory solves a similar problem. For 
instance, the form field kinetic term absorbs some four fermion couplings in 
the expression (T> A F) 2 . We can indeed view the fermionic additions in T> A as 
a covariantization of the flat derivative with respect to the form field A. There 
is, therefore, reason to expect that at least in theories that are known to have 
a holographic dual but which contain four fermion interactions the additional 
symmetries, such as local supersymmetry or the SL(2, Z) invariance in type 
IIB SUGRA, might arrange the fermion structure in such a way that the 
eventual higher fermion derivatives would still be manageable. 

6 First Class Constraints as Ward Identities 

In addition to the constraint 7i± treated above, we have to solve also the rest 
of the first class constraints Q^, Tii and J7"— . Generally speaking, they will 
impose gauge symmetry, diffeomorphism invariance and local Lorentz sym- 
metry on the boundary. Why this is not straight forwardly so is because, in 
the bulk, these constraints generate symmetry transformations through Dirac 
brackets. Due to the ansatz @ in Hamilton-Jacobi formalism, they will act 
differently on F, thus resulting in conditions that do not impose necessarily 
the full off-shell symmetries of the bulk theory. As far as the Lorentz and the 



16 



gauge symmetries are concerned the geometrically expected actions are ob- 
tained. In case the of diffeomorphism invariance some additional constraints 
arise. 

For instance, the fact that enters the generating functional F only 
through its field strength is sufficient to guarantee Q ^ = 0. This simply re- 
flects the fact that the boundary theory must have the same gauge symmetry 
as the bulk theory. A similar situation prevails as far as J7"— is concerned: 
UaJ— = just because n^L, - = 0, and the rest of the components generate 
the expected action on vielbeine and spinors, through 

V k = L V aJ ^-U±r»<p + . (55) 

This constraint guarantees, therefore, local Lorentz invariance on the bound- 
ary. Choosing gauge and Lorentz invariant S<a v will be enough to avoid 
anomalies in the boundary theory. 

The situation is somewhat more involved when the constraints that guar- 
antee diffeomorphism invariance are considered. This would mean that the 
effective action be invariant under translations generated by a vector field x- 
Or, in other words, that shifting the fields A and (p infinitesimally by their 
Lie-derivatives 

C X A = L x F + di x A (56) 

^ = (Ac + \ (57) 

we get a contribution that combines together with a contribution from the 
integration measure to a total derivative of the Lagrangian. Note that a 
general spinorial Lie-derivative does not obey the Leibniz ruleQ It is there- 
fore useful to restrict to Killing fields ^uXj) = 0, for which the formula (57) 
applies. 

Solving the constraint Hi = we get, again, assuming p odd and the 
Weyl decomposition, that the variation 

= v «' i ^+4 ¥i|! '"^ fi * +i) '4) (58) 

should annihilate the effective action. This differs from a Lie-derivative in 
two respects: First, the transformation of the form field is accompanied by 
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For a general introduction to spinors and geometry see for instance 
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a gauge transformation. Second, its action on fermions is modified by 



A x p = -t x F A Y A M-\ V (59) 
A x if = c^CpT^M- 1 (60) 

where A x = 5 X — C x . If we want to restrict to theories where the boundary 
diffeomorphism invariance still prevails, we have to put this difference to zero. 
A solution of A x Sa v = is ensured imposing the following conditions: 

1) The Clifford action of any differential form K appearing in the fermion 
couplings (pK^T A (p commute with that of i x F, i.e. 

[^xVi=0. (61) 

For instance, for = F A , F A being of odd rank, this is clearly true. 
This condition means then, geometrically, that and F^ should be 
aligned in a certain way. In addition to this the couplings should satisfy, 
in the notation of formula (70), 

(M-\ = CM~ l C . (62) 



2) If there is a kinetic term on the boundary, such as (fptp — (fl<p)tp the 
following restrictions be true 

C X F = (63) 
i x F iA Diip = (64) 

These are strong requirements, as they concern the boundary fields and 
not only couplings, and therefore really restrict, from the bulk point 
of view, the set of acceptable initial conditions. The simplest way to 
solve them is naturally to exclude the kinetic terms from the action, 
cf. end of Section [7|. However, this might be too drastic a solution, as, 
in the above, we are assuming that there exist a Killing field x\ after all, 
we are considering an interacting theory with physical sources. More 
interestingly, these conditions can be solved by assuming i x F = 0: This 
means that the Killing isometry only changes the field A by generating 
a gauge transformation. With this understanding it is sufficient to set 

C X A w . (65) 
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This would mean that there are no restrictions on the fermion fields, 
whereas the form field potential is frozen to configurations covariant 
under flows generated by x- 

Then, considering diffeomorphism invariant terms in Sd.br, one obtains the 
Ward identity 

C X W + A x v(0^) - (O v )A x p = -A x S div . (66) 

The RHS is the failure of the counter terms to satisfy the constraint Tii = 0, 
while the vacuum expectation values would signal the spontaneous symmetry 
breaking of Lorentz symmetry in the dual QFT. This is analogous to the 
analysis in Sec. ^ where scalar fields coupled to gravity were considered. 



7 Local Expansion in the Bracket Equation 

Let us consider a particular solution of the classical equations of motion that 
behaves at the boundary t — > as 

9ij (x,t) = r 2 g %{x)+0{r l ) (67) 
A A (x,t) = t p A\{x) + 0{t u+1 ) (68) 
i>(x,t) = t a ip R (x) + 0(t a+1 ) (69) 

in the coordinate system of Eq. ([TT|). The quantities with superscript R 
refer to expressions that have a finite and nonzero limit at t — ► 0. We will 
actually not need to show whether the full coupled system has solutions with 
this particular asymptotic behaviour. This would also be quite difficult — it 
was found, for instance, in [17, |3H| that free fermions scale as a = d/2 — m, 
where m is the bulk mass. In our case the coupling of the form field Fa to 
the fermions gives rise to an effective mass term. We cannot, however, fix 
the field Fa in any useful way in order to analyse conclusively the scaling 
behaviour in this coupled system. 

Instead, we impose the above scaling behaviour for some set of critical 
exponents v and a and then derive from this - and the assumption that the 
holographic dual exist at all - consistency conditions for both the bulk and 
the boundary theories. Let us consider, in particular, the case v = —n + 1 
and o = 1/2. This assignment of critical exponents has the virtue that 
the expansion of S^iv will look like an expansion in terms of the naive mass 
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dimension relevant to supersymmetric models. This choice will turn out to 
be a convenient book-keeping device, but the results we eventually get are 
valid more generally. In order to solve the full bracket equation, we will 
have to arrange the coefficients of terms that not only have the same scaling 
behaviour, but also the same structure, to cancel. So, in principle, one can 
check a posteriori for which range of scaling exponents the terms we neglected 
are still subleading and our results continue to be then valid. We shall further 
assume that all couplings such as M and Z are marginal operators, and as 
such t-independent. This assumption is not very restrictive, not even in the 
presence of scalar fields, that would render the couplings dynamical. We can 
now write a local ansatz for that is of the same functional form as (f46|) 

Sdiv = Jd d x^/§ (hR - A + h 2 F A F A + F A <p(r A <p 

+ ^Mflv-^(pv)Mv + vL A T A ^ ■ (70) 

As we shall later restrict to a scalar coupling L, it turns out that no four 
fermion terms are needed. The couplings fcijA^C?^ an d L A T A need not 
be constant on the boundary, but they are assumed marginal, i.e. time- 
independent. 

At leading order Sdiv will diverge as t~ d when t — > 0. We shall consider 
the three lowest order contributions to the equation H± = 0, or, (Sdi v , Sdiv) = 
C. 

1) At the leading order we get 

The boundary cosmological constant is therefore essentially the square 
root of the bulk one and sets the scale. We notice that depending on the 
sign of the bulk gravitational constant we can choose to look at either 
time-like or space-like boundary surfaces. However, only the space-like 
surfaces admit an asymptotically anti de Sitter solution. 

2) The 0(t~ d+1 ) equations depend on the details of Z. The equations 
simplify assuming that Z is odd in Clifford matrices and that L is a 
scalar, i.e. a flavour matrix, as we then get 

= Z + r)LM- x ZM~ x L . (72) 
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If Z is even, we get 

ZM~ l L + LM~ X Z = 0. (73) 

In particular, we could consistently set L = 0. This would mean that, 
for instance, a bulk mass term does not automatically lead to a mass 
term on the boundary. 

3) At the next order 0(t~ d+2 ) we get from the Einstein-Hilbert and the 
form field kinetic terms conditions for the couplings k\ and ki 

_±_ K ^_ rtfe _J_ = . (75) 

For d > 2 equation (|74|) allows one to compute the coefficient of the 
R term, as expected. In d = 2 this equation does not arise as \^gR is 
marginal, and it is not to be included in S&y. Instead, the impossibility 
of canceling it in the descent equations translates to the Weyl anomaly 
in d = 2 as in Eq. ([19]). 

Equation (|75|) gives the value of &2 for d ^ 2p. Note, however, that 
for middle- dimensional form fields d = 2p this equation will not have 
solutions. In the case of free form fields, as they are always marginal 



2q|, it is just the anolgue of ([74]) for pure gravity: In this case F is 
a marginal operator and it contributes to the matter part of the Weyl 
anomaly. In the interacting case the treatment of this term depends on 
the scaling dimension v. If v < then the contribution should really 
be included in the divergent part, and a middle-dimensional form field 
would not allow consistent solutions for the bracket equation. 

4) Assuming that M and L are constants on the boundary, the fermion 
terms yield the constraints 



11 K 2 kM = M + rjLM^L (76) 
k 2 A( = C + vLM^CM^L . (77) 



2d- 1 
V 



2d -2 

In all of the equations ([H]) and QT4] ) - (]77| ) we see that the relationship 
between the bulk and the boundary couplings is essentially a scale factor k 2 A. 
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It is interesting to note that the role of L is to mix bulk fermion flavours into 
new combinations on the boundary. Assuming L proportional to M would 
lead to scaling the fermion field. Furthermore, setting L = ^/—t]M the 
dynamical part to the fermion action drops out completely M = ( = 0, and 
the only contribution comes from the L-term itself 



As discussed in the previous section this solution would not produce a dif- 
feomorphism anomaly, even without restricting to symmetric form fields, 



We have investigated the relationship between diffeomorphism invariant the- 
ories and their holographic duals, showing in particular that, in a theory that 
contains fermions, nontrivial consistency conditions arise. These conditions 
restrict, for instance, the couplings of even rank form field field strengths to 
fermions. 

After explaining how the Hamiltonian reduction is performed, with par- 
ticular attention to fermion fields, in Sections || and ^| we have discussed 
how the holographic Callan-Symanzik equation and other Ward identities 
following from first class constraints get modified in the presence of fermions 
and forms. Although the gauge and the Lorentz constraints did not lead to 
surprises, the Poincare constraints resulted in an anomalous contribution in 
the diffeomorphism Ward identity on the boundary. However, one can get 
rid of these terms, either imposing conditions on the sources, as for instance 
the boundary gauge potential to be constant in the sense of Eq. (pop, or by 
the choice of a suitable counter term (|78|) . 

We have also derived relationships between the bulk and boundary cou- 
plings, finding that the role of the cosmological constant on the boundary is 
to set the scale of all boundary couplings w.r.t bulk couplings, as expected. 
Moreover, we were able to find the first terms of the expansion of the counter 
term action Sdw, with couplings fixed in terms of those appearing in the bulk 
action. This analysis shows also that one can avoid including fermion kinetic 
terms in S^w that would have a nonzero anomalous diffeomorphism variation. 




(78) 



cf. Eq. (^). 



8 Conclusions 
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There still remain open problems. One should clarify, for instance, what 
kind of restrictions higher fermion couplings impose on the duals, or what 
role other bulk symmetries, such as supersymmetry, might play in the bracket 
structure. In particular, it would be interesting to include dynamical scalar 
fields and Rarita-Schwinger fermions in these considerations. One would also 
like to know how robust our results actually are when the scaling behaviour 
of the bulk fields at the boundary are varied. Finally, the formal structures 
arising here are quite intriguing: It would be interesting to find out whether 
the (•, -)-bracket or the P-operators have a geometrical meaning in solving 
Ward identities. 
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Greek indices /x, u, . . . refer to the coordinate directions in the bulk, under- 
lined Greek indices a, (3, . . . are flat indices in the bulk, lower case Latin 
indices refer to coordinate directions in the boundary, upper case 

Latin indices A,B, . . . refer to normalized multi-indices in the bulk, and hat- 
ted upper case Latin indices A, B, . . . refer to normalized multi-indices in the 
boundary. If there is danger of confusion, symbols with a tilde, such as V, 
are used to refer to bulk quantities and symbols with a hat, such as p, to 
boundary quantities. 

We can express the connection between an arbitrary flat coordinate basis 
{e„} and a basis {e ,ej} that involves the direction of the evolution coordi- 
nate (bulk direction) eo in terms of the vielbeine 



A Notation and useful formulae 



e ~ = Nn- + N i L i 



a 



(79) 
(80) 
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Due to the algebraic constraint L i ~n a = we have e, • n = 0. For other 



properties of this frame see for instance [45|. The boundary metric and 
vielbeine are related through 

9ij = L i -L j -T]ap (81) 

r)ap = L ig V p + rjUanp , (82) 

and the boundary gamma matrices are defined by 

T n = nj:^ (83) 

T = £* a r a . (84) 

Given the Levi-Civita connection V in the bulk we can construct a metric 
connection on the boundary by setting 



V X Y = V X Y + rjn(Y ■ V x n) - rjS7 x n{Y ■ n) (85) 

for arbitrary vector fields X, Y. This connection enjoys the properties 

Vin = (86) 
n ■ ViCj = . (87) 

The spin connection in the bulk can be expressed in terms of that on the 
boundary using 

^iaf3 = Tjifc Lp — l!y diL k a — 1] Up diUa 

+77 K {j (ria 4 - n t L D (88) 
woo^ = (djN + r/N* Kij) (n« L?p - rip L„) 

-VyiV fc ] L^L'p + rj d t np] + L i[s d t L 3 ^ . (89) 

The extrinsic curvature is 

Kii = (&9<i ~ V«JVj - VjNi) . (90) 

From the point of view of the Lagrangian formalism, the momenta are 
just notation for expressions involving fields and their derivatives 



= 2n i3 Lf --T] V§ 4>M{T n , T^}tfj Up (91) 
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7T« = 


_v_ 

K 2 


(g lj tiK-K ij ) 


E A = 






X = 


1 




X = 


1 





(92) 
(93) 

(94) 
(95) 

The Poincare constraints consist of three parts W = Tij + Tljj + Ti.j H . In 
the pure gravity sector we have 

nf = ^(_J ifl+ A)-5(t r tf- I l T ( tr nr) (96 , 

«5 = -Vi(P ia 4), (97) 
where the gravitational momenta have been shifted according to 

IF = ir ij - ^g ij G (98) 

pja = pja _ L ia G _ (gg) 

In the form field sector 

nh = VI ( ~7^F A F A - F A J A 



2A 2 A 

( E A - V9J° A ) {E A ~ VtJ° A ) (100) 
Hi, = F^(E A -^J oA ), (101) 



the different signs in front of the fermionic dynamical ( |101| ) and background 
( |103| ) currents is not a surprise, as a contribution of the first one has been 



used in the definition of the electric field E A . 
In the fermion sector 



Hjjj = \^g[-mH + m)M^-2i,Z A YU) (102) 

n\ n = (i> v n Mb 1 ^ - (L>v)r n MV> + 2^z i A v n v A ^j (103) 
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The action of the covariant derivative on spinors is, by definition 

= r^ + iw^rs^ (104) 

(^) = (a^-^^r^)^, (105) 

so that 

V,(xT^) = (px)^ + Xp*P ■ (106) 

In addition to the Poincare constraints there are also constraints that 
generate the gauge transformations A — > A + dB and gauge transformations 
on the frame bundle (local Lorentz transformations) 

g A = diE iA (107) 

jap = p i\a L £\ + _ v ^ ^M{r n , T^}lfj . (108) 



B Clifford algebra 

The formulae in this appendix are taken mostly from Ref. [47]. 
B.l Even dimensions 

Consider a metric with signature r\ ah = {(— ) d - , (+) d+ }, such that d = d_ + 
d + = 2m. The Clifford algebra is span by 

{T\T j } = 2g ij . (109) 

All representations are unitarily equivalent, and the intertwining operators 
are 

rj = ATiA- 1 (110) 

_rf = c- l TiC (in) 

-r* = D- x YiD (112) 
r* = D^ViD (113) 

where D = CA T and D = TCA T . The chirality operator is 

r = r 1 • . . . ■ v d (ii4) 
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We have the following phases 

A = aA^ (115) 

C = f)C T (116) 

DD* = 5 (117) 

DD* = 5 (118) 

where \a\ = 1, f] = ±1, and 5* = 5. The Dirac conjugate is defined as 
ip = ip^A, and the charge conjugate as ip c = CA T tp*, or ip c = TCA T ifj*. The 
chirality matrix satisfies 

rt = (-) m ATA' 1 (119) 
r T = (-) m c- l vc . (120) 

B.2 Odd dimensions 

In this appendix we build a representation of an odd-dimensional D — d + 1 
Clifford algebra with (r n ) 2 = rj starting from a given even dimensional d = 
2m Clifford algebra. We actually only need to construct the correct Clifford 
matrix r n 

T n = y/(-) m + d -rig T . (121) 

In odd dimensions not all representations are equivalent. Instead, we only 
have the intertwining operators 

= (-) d - V AT^A- 1 (122) 
Tl = (-) m C~ l rj5 (123) 
= (-) m+d - v V- l YJ) (124) 

where T> = CA T , not to be confused with (0). There are two inequivalent 
conjugacy classes: which representations belong to which depends on m, d_ 
and rj. 

We can represent the odd-dimensional intertwinors in terms of their even 
dimensional counter parts as 



+1, 


A = 


A 


-1, 


A = 


AY n 


+1, 


C = 


T n C 


-1, 


C = 


C 



y' / ! : , ■,, (125) 

(126) 
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We sometimes abbreviate the sign (— ) d -i] by e. The operators T> can be 
expressed in terms of boundary operators 



V 


= 1 = -1 


(-) d -T) = 1 
{-) d - v = -1 


D D 

t]D -b 



(127) 



where D = CA T and D = T n D. 



C Fermion phase space 

Fermions can be decomposed in bulk and boundary components in essen- 
tially two Lorentz invariant ways, namely by using chirality or reality condi- 
tions. Here we consider only chirality conditions, with which we refer to the 
eigenvalues ±1 of y^fjT n . This is essentially the only Clifford matrix whose 
eigenvalues we can consider without breaking Lorentz invariance explicitly. 

We have to divide our analysis in two cases depending on the sign rj and 
details of the bulk metric. Defining ip± = ±ip± we get 

ij ± ^jT n = ±87] ij ±} (128) 

where e is the sign that appears in the relation of the Clifford matrices to 
their Hermitian conjugates ( pL22| ) . This means that the the Dirac dual of 
a spinor of definite chirality is either of the same or the opposite chirality; 
as a consequence, the Lagrangian fields that correspond to the phase space 
coordinates will be different. It is useful to note that 

en={ { ~ )d ~ d+1 ° dd (129) 
| r] d + 1 even 

Case I. Assume erj = 1. The kinetic term separates into 

-y/r^i^-Mip^ + i+Mip + ) - d t G (130) 
-^_$(V^M) V- - \i>+d t {^gM) V> + (131) 
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where 



(132) 



The fermionic phase space consists therefore of the symplectic pairs (tp, n) 
and (tp, 7r), where (p = tp- and <p> = and 



7T 
7T 



The last two terms in ( |131| ) produce a term 



^9 ij G d t9ij 



(133) 
(134) 



(135) 



in the action, and therefore cause a shift in the gravitational momentum. 
Case II. Assume erj = —1. The kinetic term separates into 



\,'^i)(v Mc. + i)+Mi)_) - i),G 



(136) 
(137) 



where 



G 



(138) 



The configuration space is span by <p = ip_ and (p = and the momenta 
are 



7T 
7T 



(139) 
(140) 



Notice that, due to ( |125| ), the Dirac conjugates <p and ip- are formed differ- 
ently: the former using the matrix A and the latter with A. This will result 
in an extra T n everywhere, including the kinetic term, and the resulting extra 
sign hence cancels out. The gravitational momenta are shifted as in Case I. 
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